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The  concept  of  map  costyar ison  is  treated  in  a  yeneiai  ^uu 
foL  thematic,  regularly  yridaed  naps.  Three  methods  of  nay 
comparison,  direct  methods,  trend  surface,  and  £(atial  auto¬ 
correlation,  are  investigated.  ibe  direct  methods  and  treru 
surface  method  are  applied  tc  artificial  surfaces  anu  sultt- 
guently  to  regularly  gridded  maps  of  the  pollutant  total 
suspended  particulates.  A  method  of  forming  spatial  :r  twe 
dimensional  a u tocor rela tiou  maps  is  presented.  A  strategy 
based  on  these  methods,  is  described  whereby  a  laiye  nuster 
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Chapter  I 
IITBODUCTIOH 


1.  1  POBPOSE 

The  purpose  of  this  investigation  was  to  compare  some  ex¬ 
isting  Methods  of  aap  coaparison  and  investigate  potentially 
useful  computer-assisted  methods  of  nap  comparison.  Ihe  em¬ 
phasis  has  been  on  comparing  the  different  methods  ana  de¬ 
veloping  a  possible  strategy  for  comparing  large  numbers  of 
potentially  similar  maps. 

There  are  many  reasons  for  comparing  maps  and  an  increas¬ 
ing  number  of  gridded  data  map  sources  are  becoming  availa¬ 
ble.  Geographers  have  compared  maps  of  many  types  of  human 
variables  such  as  heat  source  vs  income,  spead  of  infection 
vs  well  water  source,  and  many  other  sets  of  variables. 
Geologists  compare  maps  of  different  subsurface  rock  units 
to  study  tectonic  history  of  an  area.  In  pollution  studies, 
maps  of  different  variables  such  as  terrain,  wind  patterns, 
and  emissions  sources  are  compared.  For  this  in vestigation 
monthly  pollution  maps  will  be  compared  in  order  to  demons¬ 
trate  possible  computer-assisted  methods  of  map  comparison. 


t 
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The  analytical  skills  of  a  huiao  investigator  will  n 
replaced  by  a  clever  program,  but  helping  the  investi 
focus  his  attention  on  potentially  important  relation 
between  naps,  without  having  to  look  at  every  possible 
bination  of  naps.  Bakes  the  developaent  of  computer-ass 
nethods  of  aap  comparison  worthwhile.  If  the  geogra 
cartographer,  or  geologist  can  eliainate  or  reduce  th 
ount  of  initial  coaparisons  that  require  manual  eft 
then  the  investigator  is  free  to  concentrate  on  the  p 
tially  iaportant  aap  relationships. 

1.2  METHODS  OF  HAP  COHPABISOM 

In  this  investigation,  three  methods  of  sap  coapa 
are  discussed  and  demonstrated.  The  naps  that  are  cca 
are  thematic,  regularly  gridded  arrays.  Computet  pro 
were  written  to  implement  the  comparrson  metnois.  In 
testing  and  evaluation  of  the  programs  was  done  witn  at 
cally  generated  maps.  A  final  demonstration  ot  the  com 
son  methods  was  done  using  pollution  data  obtained  fro 
Ohio  EPA  [ Ohio  EPA,1979J. 

Hap  comparison  methods  used  in  this  investigation 
into  three  general  catagories:  direct  methods,  trend 
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ods,  tread  surface  and  autocorrelation.  Ibe  direct  methods 
coapare  the  original  nap,  or  array  sealers,  ia  a  one  to  one 
relation  aith  other  aaps  .  The  direct  aethods  include  visu¬ 
al  coaparison,  isopach  aaps,  and  direct  correlation.  The 
trend  surface  aap  coaparison  techaigues  fit  a  Li-variate  po- 
lynoaial  to  the  gridded  data  sets,  using  the  aethoa  of  least 
squares.  The  estiaated  paraaeters  of  the  polynomial  fit  are 
coapared  using  the  correlation  cofficient.  The  third  method 
uses  the  autocorrelation  of  each  aap  to  generate  a  represen¬ 
tation  of  the  surface.  The  autocorrelation  of  a  series  or 
array  is  coaputed  by  coaparing  the  series  or  array  to  it¬ 
self,  the  exact  aethod  that  is  used  to  do  this  will  be  dis¬ 
cussed  in  chapter  five-  Autocorrelation  can  be  thought  of 
as  the  internal  relationship  of  a  set  of  data.  A  limited 
autocorrelation  function  is  used  to  produce  a  autocorrela¬ 
tion  curve  for  each  of  the  aaps.  A  spatial  autocorrelation 
function  is  a  aeasure  of  the  internal  consistency  of  a  sur¬ 
face  (aap).  The  data  set  has  two  diaensions,  thus  the  re¬ 
sult  of  calculating  the  spatial  autocorrelation  function  is 
an  array.  A  spatial  autocorrelation  function  can  be  usea 
to  fora  an  autocorrelation  aap  for  each  of  the  original  sur¬ 


faces. 


The  direct  aethods  of  comparison  include  the  visual  com- 
parison  of  two  aaps.  In  order  to  visually  coapare  the  arid- 
ded  data  arrays,  one  oust  represent  the  data  in  a  aore  human 
compatible  fora.  The  Geodetic  Science  Plotting  Package 
(GSPP) ,  [ SUNKEL, 1980  ],  is  a  plotting  package  that  will  pro¬ 
duce  isarithaic  aaps  frou  regular  gridded  arrays.  GSPP  uses 
a  bi-cubic  spline  method  of  interpolation  to  torn  isarithaic 
naps  froa  regularly  gridded  arrays  of  data.  All  isarithaic 
aaps  presented  in  this  thesis  were  formed  using  GSPP. 

1.3  HAPS  OSBD  FOB  COHPABISOB 

For  the  purposes  of  this  investigation,  mappeu  da*a  is 
regarded  as  a  set  of  observations  ordered  with  respect  to 
two  spatial  coordinates.  The  spatial  coordinates  are  repre¬ 
sented  with  a  regular  grid  structure.  The  distinction  tet- 
ween  irregular  and  regular  gridded  data  can  have  rather  sig¬ 
nificant  nuaerical  consequences.  Only  regularly  gridded 
data  was  used.  Xo  avoid  the  problems  inherent  in  the  gria- 
ding  of  irregularly  distributed  data,  all  discussions  with 
respect  to  aap  coaparison  methods  will  assume  that  observed 
values  come  froa  the  intersections  of  a  regular  grid.  The 
pollution  data  used  caae  froa  irregularly  located  collection 
sites.  Ihe  aethod  and  theory  used  to  grid  the  data  will  Le 


discussed  later 


The  gridded  data  sap  has  been  referred  to  as  a  "numerical 


■ap"  or  a  "digital  map"[  EtOBINSON  et  al.;  1978J.  It  is  not  a 
particularly  informative  map  if  it  is  viewed  in  its  most 
elemental  fora.  It  would  appear  as  only  a  grid  of  points 
with  real  number  values  printed  near  the  points.  Normally, 
a  gridded  data  set  would  be  presented  as  an  isoline  map, 
that  would  have  been  formed  by  interpolating  between  the 
gridded  data  points  to  represent  isolines. 

(lore  specifically,  the  isarithaic  map  is  a  statistical 
surface.  That  is,  it  may  represent  the  volume  or  surface  or 
data  in  a  mapped  area.  The  isarithaic  map  is  of  two  general 
types  the  isoline  and  isopleth  map.  The  transition  between 
the  isolines,  of  the  isarithaic  map,  is  assumed  to  ce  conti¬ 
nuous.  On  the  other  hand,  the  isopleth  map  represents  aata 
that  can  not  exist  at  a  point.  Ihe  statistic  represented  by 
the  isopleth  map  occurs  inside  a  boundary  of  finite  extent. 
The  boundaries  often  used  in  practice  are  enumeration  dis- 
ticts,  townships,  counties,  etc.  Possible  statistics  used 
are;  average  sale  price  of  homes,  cattle  per  acre,  deaths 
per  thousand,  etc.  The  isoline  and  isopleth  maps  are  very 
different  with  respect  to  the  nature  of  the  data  they  are 
used  to  represented. 
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The  appearence  of  the  isoline  and  isopletb  naps  art  *uit e 
sinilar.  They  both  use  lines  as  the  napping  syatol.  lhe 
isopleth  nap  appears  very  sinilar  to  the  isoline  nap.  They 
both  represent  a  statistical  surface  that  is  assumed  to  Le 
continuous,  yet  the  statistic  of  the  isopleth  map  can't  ex¬ 
ist  at  a  point  because  isopleths  represent  data  that  is  av¬ 
eraged  over  an  area  [  fiOBIM-SON;  t978#  p.  225  ], 

The  nap  data  used  in  this  in vestigation  is  oi  the  type 
that  an  isarithnic  nap  would  be  constructed  from.  The  iso¬ 
lines  are  "threaded"  anong  the  gridded  data  points  ty  an  ap¬ 
propriate  algorithm.  The  result  is  an  isanthmic  nap  that 
is  a  visual  representation  of  the  gridded  data.  It  would  it 
well  to  enphasize  that  the  techniques  used  in  this  investi¬ 
gation  assune  the  data  to  be  in  the  form  of  regjlai  guddtu 
values,  not  stored  as  isarithns.  It  is  possible  to  overlook 
the  distiction  between  an  isarithnic  nap  that  was  formed 
fron  gridded  data  and  the  gridded  data  set  it  was  loused 
from.  The  gridded  data  sets  are  what  Mill  be  compared. 

Because  this  investigation  of  nap  comparison  methods  is 
linited  to  regularly  gridded  data,  it  would  be  a  fair  ques¬ 
tion  to  ask,  "where  would  one  obtain  data  in  this  special 
fornat  and  what  value  are  the  techniques  that  address  only 


regular  gadded  data?".  To  answer  the  first  part  of  thrs 
question,  there  are  a  number  of  sources  of  qridded  data. 
The  digitizing  of  existing  aaps  has  produced  very  large  data 
bases  of  regularly  gridded  elevation  data.  These  data  bases 
are  often  referred  to  as  digital  terrain  aodels  (Din)  or  di¬ 
gital  elevation  aodels  (DEH).  The  Defence  Happing  Agency 
(DBA)  has  a  very  extensive  DTH  £ SAULH1EB; 1979  ].  Landsat  Im¬ 
agery  is  collected  directly  in  the  gridded  (raster)  forsat 
and  is  readly  available.  Another  potentinally  important 
source  for  regularly  gridded  data  is  the  analytical  phcto- 
graaaetric  instruaent  [ HELAVA ; 1 978 J.  A  lesser  source,  but 
still  an  iaportant  one,  is  the  collection  of  geophysical 
data.  Electrical,  gravity,  and  nagnetic  surveys  are  often 
collected  on  a  regular  or  seal-regular  grids.  Finally,  a 
potentially  large  source  of  gridded  data  is  the  transf or  na¬ 
tion  of  raadonly  or  irregularly  distributed  data  to  a  regu¬ 
lar  gridded  format. 

The  second  part  of  the  question,  "what  value  are  techni¬ 
ques  that  address  only  regular  gridded  data?"  is  the  more 
interesting  aspect  of  the  question.  There  are  a  number  of 
uses  that  have  been  described  for  this  type  of  data.  Oue  of 
the  primary  uses  is  for  navigation.  Hap  Beading  is  a  stan¬ 
dard  navigation  technique  where  the  navigator,  from  some 
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vantage  point,  compares  the  local  terrain  with  a  chart,  then 
orients  hiaself,  and  finally  determines  his  precise  loca¬ 
tion.  The  carrent  version  of  this  nethod  replaces  the  navi¬ 
gator  with  a  conputer  and  radar;  the  chart  is  replaced  by  a 
DTH.  The  terrain  is  scanned  by  the  radar  system  under  com¬ 
puter  control  and  the  DTH  is  coapared  to  the  terrain.  Until 
DHA  released  its  DTH  to  the  civilian  agencies,  research  that 
reguired  a  dense  DTH  was  difficult.  The  DTH  has  been  com¬ 
pared  to  geological  inforaation,  pollution,  and  wind  pat¬ 
terns  [TESCHE  and  BEBGSTBOR  ; 1978 J,  to  name  a  few.  Hot  only 
can  aaps  of  different  types  he  coapared,  but  gndded  data 
sets  that  vary  temporally  can  also  be  studied. 
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1.«  H1I54J.  HAP  SOMPiiisei  fll IflflflS 

flanual  aethods  of  map  coaparison  ace  cot  efficient  enough 
to  handle  lacge  aaounts  of  data  oc  Baps.  In  order  to  com¬ 
pare  a  large  nuaber  of  surfaces,  the  investigator  Bust  have 
soae  aethod  of  reviewing  the  relationship  between  all  possi¬ 
ble  aap  combinations.  The  visual  coaparison  of  maps  is  a 
useful  tool,  but  the  nuaber  of  coaparisons  for  n  traps  is: 

t  =  0.5  (n*2)  (n*  1)  , 

t  :  total  possible  aap  pairings 
n  :  nuaber  of  aaps  to  be  coapared 

when  all  possible  coabinations  are  considered.  As  a  prelim¬ 
inary  coaparison  aethod,  the  visual  aethod  does  not  appear 
to  be  very  useful.  The  other  direct  aethod,  isopach  maps, 
is  also  inefficent,  because  it  produces  so  many  difference 
aaps  that  it  is  difficult  to  understand  what  the  relation¬ 
ships  are  between  the  original  aaps. 

A  coaputer-assisted  aethod  of  coaparison  seems  to  re  ne¬ 
cessary  for  preliainary  stages  of  investigation.  Two  compu¬ 
ter-assisted  aethods  of  preliainary  coaparison  are  used  in 
this  investigation,  they  are  direct  coaparison  of  tbe  gnu- 
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ded  data  sets  aad  comparison  of  parameters  of  trend  surfaces 
of  data  sets.  These  methods  are  applied  to  artifical  sur¬ 
faces.  The  relationship  between  the  artifical  surfaces  was 
known.  Once  the  known  surfaces  were  compared  and  a  "base¬ 
line"  relationship  between  the  artifical  surfaces  was 
formed,  the  next  step  was  to  apply  the  cooputet-assisted 
methods  to  ata  that  had  an  unknown  relationship. 

The  autocorrelation  of  a  surface  say  also  have  potential 
as  a  nap  comparison  tool.  The  autocorrelation  function  of  a 
surface  is  unique  and  could  be  used  to  represent  a  map.  1  he 
comparison  between  the  autocorrelation  function  of  two  naps 
could  give  an  insight  into  the  relation  between  the  gridded 
data  maps.  Obviously,  autocorrelation  will  net  re  practical 
as  a  preliminary  comparison  method,  but  it  should  it  usetul 
as  a  secondary  method  of  sap  comparison. 

Hhile  the  methods  of  comparison  used  are  not  a  complete 
list  of  every  possible  technique,  they  do  include  a  wide 
range  of  theory  and  application.  The  trend  suirace  method 
has  been  used  in  the  past  but  with  much  smoother  surfaces. 
To  the  best  of  the  authors  knowledge  the  spatial  autocorre¬ 
lation  function  has  not  been  used  before,  as  a  sap  compari¬ 
son  technique,  comparison  technigue.  For  the  map  comparr- 
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son  aethods  that  ace  to  be  demonstrated  the  regularly 
gridded  froaat  is  the  aost  pracitical. 

Chapters  three,  four,  and  five  describe  the  nap  compari¬ 
son  aethods.  Chapter  six  describes  the  artifical  surfaces 
that  were  used  to  evaluate  the  computer-assisted  methods  of 
aap  comparison.  Additionally,  chapter  five  gives  a  descrip¬ 
tion  of  the  pollution  data  aaps  that  are  to  oe  compared. 
Chapter  seven  presents  the  results  of  the  nap  comparison 
aethods  applied  to  the  artifical  sufaces  and  the  pollution 
aaps.  Chapter  eight  contains  a  strategy  for  coiputei-as- 
sisted  aap  coaparison  using  the  methods  that  have  been  dis¬ 
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Chapter  II 

8ETIEH  OP  BBS.A1E0  B8SEABCH 

Hap  comparison  techniques  are  used  whenever  some  type  of 
spatial  phenomenon  is  to  be  investigated.  Any  discipline 
that  uses  spatial  data  has  a  need  for  methods  of  comparing 
maps,  with  the  fields  of  geology  and  geography  notable  in 
this  respect.  Naturally  earth  scientists  have  seen  leaders 
in  the  area.  Hany  of  the  techniques  used  in  this  paper  are 
drawn  from  these  disciplines. 

An  early  compter-assisted  method  of  map  comparison  was 
described  by  fiobinson  and  Bryson  £1957].  They  us'd  orthogo¬ 
nal  polynomials  to  calculate  best  fit  east  west  pronles  of 
population  density  and  average  rainfall  for  the  state  of  Ne¬ 
braska.  The  parameters  of  the  univariate  pol ynomials  were 
then  standarized  and  a  correlation  catticieut  was  calcualt- 
ed.  The  study  compared  a  profile  from  a  isopleth  sap  with  a 
profile  from  an  isarithmic  map.  The  resultant  correlation 
cofficient  comfrimed  the  "graphic  correlation"  that  had  been 
observed  between  population  and  average  rainfall.  The  next 
advance  in  this  technique  cane  from  the  field  of  geology. 
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The  comparison  and  contrasting  of  geological  features 
from  one  area  to  another  is  a  fundamental  interpretive  skill 
needed  by  all  geologists.  A  somewhat  narrower  version  of 
the  same  skill  is  the  comparison  of  maps  iron  the  sane  area 
with  different  z-values.  More  specificity,  a  geologist  may 
wish  to  compare  maps  of  "different  chemical  or  mncr ulo^ical 
constituents,  different  grain  size  parameters,  or  structual 
maps  on  different  geological  horizions"  [DAVIS,  1973,p.37uj. 

An  example  of  the  comparison  of  structural  maps  is  the 
now  classic  paper  by  Herriaa  and  Sneath  [196b].  A  structur¬ 
al  geological  nap  relates  the  depths,  to  the  top  or  a  hori¬ 
zon,  with  some  datum,  usually  sea  level.  Drill  ncl?  data  is 
usually  the  source  of  the  depth  information.  For  the  Sctn- 
an  and  Sneath  investigation,  the  structual  maps  used  were  of 
subsurface  beds  that  extended  throughtout  Kansas.  Iht. 
structual  surfaces  that  were  compared  were  treated  as  mais 
to  be  compared.  Ihe  reason  for  the  study  was  to  show  that 
some  groups  of  beds  had  experienced  similar  tectonic  forces. 
That  is,  the  characteristics  of  the  surfaces  should  be  simi¬ 
lar  had  they  undergone  deformation  at  the  same  time. 

The  method  of  map  comparison  used,  involved  the  fitting 
of  low  order  polynomials,  by  least  sguares,  to  randomly  se- 
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lected  sites  on  structual  geological  aaps.  Estimated  param¬ 
eters  of  the  polynomials  for  each  surface  were  usea  to  cal¬ 
culate  a  correlation  coefficient.  the  correlation  coefti- 
cient  coaputed  in  this  aanner  represents  the  correlation 
between  respective  normalized  parameter  values.  Ej  calcu¬ 
lating  the  correlation  coefficients  for  each  possible  pair¬ 
ing  of  surfaces,  flerriaa  and  Sneath  were  tetter  able  to  un¬ 
derstand  the  similarity  between  the  maps. 

Both  geology  and  geography  use  the  above  technique  of 
comparing  estimated  parameter  values  of  different  trend  sur¬ 
faces  [CLIFF,  1975].  Another  map  comparison  technique  that 
is  shared  by  both  disciplines  is  the  isopach  map  method. 
With  this  method,  the  maps  that  are  to  be  compare!  ate  at¬ 
tracted  from  one  another.  The  result  is  a  third  map  mat 
represents  the  difference  between  the  original  surfaces. 
This  is  a  very  common  technique  [DAVIS,  1973;  NChniCK, 
1972].  A  major  shortcoming  of  this  method  is  the  neon  to 
interpret  the  isopach  map.  That  is,  the  resulting  isopach 
maps  don't  give  an  easily  quantifiable  measure  of  the  simi¬ 
larity  between  the  two  initial  maps. 

Additionally,  there  are  a  set  of  techniques  that  are 
found  in  geology,  geography,  and  geodesy  that  are  quite  si- 


■liar  in  nature  but  very  different  in  their  applications. 
These  techniques  are  kriging,  autocorrelation ,  and  colloca¬ 
tion.  They  all  are  used  to  study  the  spatial  relationship 
of  soie  variable.  In  geography,  autocorrelation  could  ne 
used  to  study  the  spread  of  aeasles  infection^ CLIP F  and  GdD, 
1981  ].  The  kriging  aethod  of  geology  is  used  to  predict  the 
spatial  distribution  of  ore  grades£ CLAEK ,  1979;  kcval  ct 
al«,  I960].  In  geodesy,  point  and  lean  gravity  anouliies 
are  the  variables  predicted  with  collocation  [SUNKF.L,  1981], 
More  specificlly,  the  correlograa  of  kriging  ana  the  covari¬ 
ance  function  of  collocation  have  a  similar  tore.  Tnese 
functions  are  in  turn  auch  like  the  autocorrelation  func¬ 
tion-  They  relate  a  saapled  quantity  at  one  location  with  a 
siailar  quantity  saapled  at  a  location  a  distance  ^a)  away. 
An  average  value  of  this  relationship  is  calculated  tor  all 
samples  that  are  located  (d)  distance  apart.  Then  a  surilar 
calculation  is  done  for  distance  (2d)  and  similarly  with  in¬ 
creasing  distances  (nd).  The  relationship  used  in  all  three 
methods  is  the  variation  between  the  saaples. 

The  spatial  autocorrelation  function  of  geology 
[AGTEBBEBG  ,1974]  seeis  to  be  the  aost  promising  method  tor 
aap  coaparision.  The  values  of  the  autocorrelation  function 
have  a  liaited  range  of  -1<r<1  where  as  the  variogram  as 
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used  ia  Ktiging  will  have  values  that  are  dependent  on  the 
variance  of  the  variable  osed.  The  theory  and  use  of  the 
autocorrelation  function  as  used  in  geography  is  priuci^iiy 
associated  vith  ordinal  and  interval  data  £  CLIFF  et  al., 
1975;  CLIFF  and  OBD#  1978],  but  the  spatial  autocorrelation 
function  as  in  geology  and  the  covariance  function  of  geode¬ 
sy  seen  to  be  the  lost  practical  functions  foi  map  compari¬ 
son.  They  are  easily  applied  «.o  gridded  data  naps  and  each 
surface  should  have  a  unique  autocorrelation  function. 


i 
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Chapter  III 

DIfiECT  H BIHODS  OF  BAP  COBPABISOB 


3- 1  USUAL  COBP ABISOI 

Possibly  the  siaplest  aethod  of  aap  coapansor.  is  to  form 
isarithaic  maps,  froa  gridded  data  sets  (the  maps  or  surfac¬ 
es  of  this  investigation)  and  coapare  them  silo  Ly  side. 
This  direct  aethod  is  often  done  as  an  initial  step  prior  to 
other  types  of  investigation  or  coaparisons.  A  direct  cou- 
parison  has  soae  practical  probleas,  but  is  worthwhile  as  a 
limited  use  technique. 

Hays  of  coaparing  two  aaps  visually  include  examining 
then  side  by  side.  Baking  transparent  overlays,  or  using  a 
light  table.  A  light  table  was  found  to  be  the  most  practi¬ 
cal  aethod  for  viewing  the  aaps  together.  For  tnis  to  worn 
well,  the  scale  of  the  two  aaps  should  be  the  same,  and  si- 
ailar  vertical  intervals  should  be  used.  The  following  are 
the  aost  coaaon  things  to  look  for  with  this  aethod. 
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1.  Positive  and  negative  similarity,  do  tbe  peaks  and 
depressions  of  the  surfaces  coincide? 

2.  Do  isolines  cross,  remain  parallel  or  at  least  nearly 
parallel? 

3.  Does  an  isoline  of  a  particular  value  on  one  map  fol- 
Iom  the  same  general  pattern  as  an  rscline  of  the 
same  value  on  the  second  map? 

□sing  the  visual  aethod  has  obvious  linitations.  Comparison 
of  aore  than  a  fee  saps  quickly  becomes  tedious.  It  is  dif¬ 
ficult  to  reaain  objective  when  comparing  many  different  map 
combinations.  The  map  that  is  on  top  has  a  tendency  to  do¬ 
minate  the  comparison.  The  author  has  forme!  difterent 
judgments,  on  the  similarity  of  two  different  maps,  after 
the  maps  were  reversed  and  viewed  a  second  time.  A  map  with 
dense  isolines  in  one  area  can  easily  draw  attention  away 
from  similarities  that  occur  elsewhere  in  the  two  maps. 
This  method  was  applied  to  the  pollution  maps  and  spatial 
autocorrelation  maps.  The  success  of  the  visual  ccmpar ison 
method  will  be  discussed  in  chapter  seven. 

The  problems  inherent  in  this  type  of  direct  comparison 
are  strong  arguments  for  more  objective  and  automated  forms 
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of  aap  comparison.  The  rest  of  this  chapter  will  discuss 
other  direct  foras  of  aap  comparison,  direct  correlation  of 
points  in  arrays  and  isopach  aaps. 

3.2  SI H F1.E  CORBBLiTIOl 

With  this  aethod,  the  correlation  between  two  sufaces  is 
developed  without  regard  for  saaple  location.  The  correla¬ 
tion  between  z-values  at  corresponding  grid  points  ,  on  two 
different  naps,  is  calculated  using: 


rJh:  correlation  cofficient  for  surfaces  j  and  k 
COVJk:  covariance  between  surfaces  j  and  k 
Sj  :  standard  deviation  of  surface  j 
n:  the  total  nuaber  of  grid  points  on  the  surface 


/  n) 


where 
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Obviously,  this  aethod  has  one  serious  drawbacks.  The  fun- 
daaental  property  of  a  aap  (spatial  location)  is  alaost  coa- 
pletely  ignored.  Id  a  very  general  way,  the  relationship 
between  two  sufaces  is  represented,  but  the  spatial  inforaa- 
tion  is  oot  considered.  The  data  points  are  treated  as  ob¬ 
servations  without  spatial  inforaation.  The  aaps  are  com¬ 
pared  as  though  they  were  a  series  of  points  Hatched  to 
another  series  of  points. 

3.3  ISOPICa  METHOD 

This  aethod  produces  a  aap  that  represents  the  differenc¬ 
es  between  two  aaps.  The  corresponding  locations  on  two 
aaps  are  subtracted  froa  each  other. 


**  =  UiK  -  \ 

Djk:  the  jk  th  eleaent  of  aap  Q. 
vJK :  the  jk  th  eleaent  of  aap  V. 
Ijk:  the  resultant  eleaent  of  the 


Isopach  Hap. 
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The  isopach  aap  that  is  constructed  will  show  isolines  of 
the  differences.  This  net hod  is  by  illustrated  using  waps 
that  are  the  saae  except  that  noise  was  added  to  one  of  the 
surfaces.  The  aethod  of  adding  noise  will  be  discussed  in 
chapter  seven.  The  resnltant  nap  is  a  graphical  representa¬ 
tion  of  the  noise  < see  figure  1). 


figure  1:  Isopach  Bap;  Noise  without  Trend 
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The  simplest  use  of  an  isopach  aap  in  aap  coapanson  is 
to  look  foe  soie  type  of  trend  in  the  nap.  That  is,  to  look 
for  structure  in  the  surface.  Zf  no  structure  is  apperant 
as  in  figure  1,  then  the  difference  could  be  assumed  to 
represent  only  random  variations.  The  priaary  advantage  ot 
this  aethod  is  that  it  is  quite  siaple  and  gives  a  guicx 
look  at  the  siailarity  of  two  saps.  The  user  should  be 
aware  that  isopach  naps  are  not  unique;  similar  isopach  maps 
aay  result  from  very  different  surfaces.  As  Davis 
has  pointed  out  an  aabiguous  cases  can  arise  with  profiles 
that  are  subtracted  froa  each  other.  Hathematical ly  correct 
differences  can  give  the  sane  surface  fron  pairs  ot  surfaces 
that  are  very  different. 

Figure  2  shows  two  dinensional  examples  of  this  prcrlec. 
These  aabiguous  cases  are  sore  likely  to  occur  when  stand¬ 
ardized  data  is  used  to  construct  the  isopach  map. 


Chapter  ZV 

TB2ID  SDBFACE  H1TH0DS  OF  BAP  COFPAIISOS 


4. 1  IBTBODflCTIOl 

Trend  surface  analysis  is  a  eethod  of  descnacg  a  set  of 
variables  that  are  related  spatially.  Using  multiple  re¬ 
gression  large,  systematic  variations  of  the  dependent  van- 
atle  are  described.  Trend  surface  sethods  are  usually  based 
on  one  of  two  aodels,  the  power  series  polynomial  or  tngo- 
noaetric  polynoaial  [ POECKEH, 1972, p.26  ].  The  General  Linear 
Hodel  is  the  unifiying  theory  for  the  two  trend  surface  ao¬ 
dels.  The  General  Linear  Bodel  has  the  form  [ GRAYEILL  and 
KHDHBEIH,  1965]  : 


k 

*  *  *>  ♦  ♦ « 

!=i 

The  X  and  1  are  observable  random  variables,  £o  and 
are  the  unknown  parameters.  The  value  e  is  an  unobservable 
or  residual  ;  it  represents  the  difference  between  the  ex¬ 
pected  and  observed  I.  It  is  the  value  that  is  to  be  mimic- 
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ized  by  the  least  squares  aethod.  The  power  series 
polynomial  aodel  when  used  in  the  general  linear  aodel  caus¬ 
es  tbe  saaaation  tern  to  expand,  giving: 

i  ♦/%  . . 

for  the  single  variable  case.  loraally  the  bi-variate  case 
is  used  for  surface  fitting.  The  general  for m  being: 

z  W/v  */v  'A1' 'A' 1  'M'  "•*AxV*  e 

where  X,  l,  and  Z  are  the  observable  variables. 

When  coaparing  naps  using  tbe  trend  surface  method  the 
estiaated  paraaeters  of  the  aatheaatical  aodel  are  the  im¬ 
portant  values.  By  computing  the  correlation  retween  the 
paraaeters  of  two  different  maps,  a  quantitative  method  of 
comparison  is  possible.  The  correlation  between  two  diffe¬ 
rent  sets  of  paraaeters  is: 

C0¥  <  B,  ,  B*  ) 
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Bj  :  the  estimated  parameters  of  surface  i. 

S|  :  the  standard  deviation  of  the  paraaeters  for  sur¬ 
face  i. 


The  correlation  coefficient  calculated  in  this  mduuei  gives 
a  useful  aethod  for  comparing  different  naps.  The  correla¬ 
tion  coefficient  has  a  known  range,  ,  which  aids  the 
coaparison.  The  following  discussion  will  briefly  describe 
the  polynoaial  mathematical  models  that  are  used,  the  method 
of  solution,  and  soae  of  the  auaerical  considerations  asso¬ 
ciated  with  fitting  a  polynoaial  to  regularly  gnddeu  data. 

4. 2  POLfOBIAI  MODELS 

The  power  series  polynoaial  used  for  treuu  surfaces  rs 
gurte  flexible  and  can  have  aany  different  forms.  The  num¬ 
ber  of  terns  or  degree  of  a  polynoaial  is  unlimited.  The 
naaing  of  the  polynoaial  models  follows  two  groupings.  lhe 
first  and  nost  coaaon  is  where  a  polynoaial  is  named  after 
its  "degree".  The  degree  refers  to  highest  power  to  which 
any  tern  is  raised.  The  pattern  is  illustrated  in  figure  3. 
The  nuaber  of  teras  (t)  in  this  type  polynomial  is 

t  =  1/2  (  a*2  )  (  n* 1  )  ,  n  :  degree  of  polynomial. 
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Figure  3:  Degree  flethod  of  Haling  (after  Kratky,l975) 

The  second  naiing  scheie  is  the  bi-polynonal.  It  is 
named  by  the  aaxiaua  power  of  individual  terns  and  not  the 
product  of  powers  £  KRATK1  ,1975].  Figure  4  shows  the  pat¬ 
tern  that  is  used  in  naiing  bi-polynomials. 
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Figure  4:  Bi-polynoiial  Baling  Scheie  (after  Kratky,l975) 
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The  teras  bilinear  (d=1  ) ,  biquadratic  (d=2  ) ,  bicubic  etc. 
are  used  to  describe  this  type  of  polynoaial.  The  nuncer  of 
teras  is  n2 . 

Three  aatheaatical  aodels  or  polynoaials  were  used  in 
this  investigation,  these  aere: 

1.  the  cubic  polynoaial  (10  teras) 

z-/?o  ♦  /s,  *  ♦  p,  /%  *3*  A  I  -  A'1' 

A*3?  ‘A,1  *  *A?!i  • 

2-  the  bi-cubic  (16  teras) 

z=  ft  o  ♦  ft,  z  ♦  A  x2  ♦  x3  ♦  /?4  y  ♦  A  >2  * 

A  y3*  Ax  7  *Ax2y  ♦  Ay2x  +Ao*2y2* 
Al3y  */V  y3  *  A*2  73  *  A*3  72  *  A5*3  y3 ' 

3.  the  Li-guartic  (25  teras) 

2=  A 0  ♦  A x  ♦  A  12  ♦  A  *3  4  A  y  +  A  y2  4 
A7 3 'A*  7  *Ax2y  *Ay  y2  +  /V2y2f 
Ax3y  4  A/  ?3>,!^Ax3y"4  AU x4* 

fipt'fi*  74 4 As*4 y4 4 A4*3 y3 * 


«.3  ifijoo  si  ssmisi 


The  least  squares  aethod  of  solution  was  used,  in  solving 
for  the  paraaeters  of  the  bi-variate  polyooaial.  Using  1a- 
trix  notation  the  general  fora  of  the  solution  is 
£  OTOLI A , 1 967 ] 


/3  *  U'P*)"1  **^1 


A: 


L: 


P: 


the  design  matrix 

vector  of  observation  for  the  dependent  z  variatle 
weight  matrix 

vector  of  parameters  of  the  polynomial  model 


Using  equal  weights  the  above  equation  reduces  to 
/?  =  (A*  A)-1  A'L. 

j  q  c  n  n  1 

n  =  m  1  o. 

'  c  c  c  1 


c;  nuaber  of  paraaeters  in  polynomial  model 
n:  nuaber  of  observations 
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The  last  equation  is  the  so  called  "normal  equations".  The 
■os t  difficult  task  in  the  solution  of  the  noraal  equations 
is  coaputinq  the  inverse  of  K.  This  inverse  can  be  rather 
troublesome  to  conpute  explicitly,  thus  the  next  section 
(4.4)  will  discuss  soie  of  the  pitfalls  associated  with  the 
trend  sarface  method. 

Two  methods  were  used  to  evaluate  the  ability  cf  the  po¬ 
lynomial  model  to  fit  the  qridded  data  sets.  The  first  was 
the  computation  of  the  <v,v>  or  sum  of  squares  of  the  resi¬ 
duals.  The  second  method  was  the  goodness  of  fit  method 
[ SERBIAN  Ah D  SNEATH;  1966  ]: 

2>?  -  <2>’!  /■ 

p  *  100 - . 

/  ° 

x,  :  the  value  of  the  trend  surface  at  location 
of  data  point 
xQ  :  observed  data  value 
n  :  number  of  data  points 


The  sun  of  squares  and  goodness  of  fit  were  calculated  with 
each  trend  surface  computed. 
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I. 4  ■aiBBIClL  C0BSIDBBATIQ1S 

There  are  three  general  sources  of  error  for  polynomial 
trend  surface  calculations:  ill-conditioning  of  the  N  ma¬ 
trix,  mathematical  error*  and  technical  error  [  Unwin,  1975  j. 
The  solution  of  the  nornal  equations 

*/3  =  u 

/?  =  H  1  D 

is  the  usual  source  of  error  in  the  trend  surface  method. 
The  parameters  aust  be  estiaated  as  accurately  as  possible 
rf  they  are  to  be  used  to  represent  the  original  data  i.hou 
comparing  the  surfaces. 

The  conditioning  of  a  aatrix  is  actually  a  measure  of  the 
sensitivity  of  the  aatrix  to  change.  If  a  snail  change  in 
the  original  aatrix  can  cause  an  unusually  large  change  in 
its  inverse  aatrix*  then  the  original  aatrix  is  said  to  Le 
ill-conditioned.  Theoretically,  the  H  aatrix  should  be  po¬ 
sitive  definite,  and  nou-singular 
[ GRAIBILL, 196 1 ; UN WIN , 1975  ].  Thus  one  can  assume  that  an  in¬ 
verse  exists,  but  to  know  boa  ill-conditioned  the  matrix  is 
soae  tests  should  be  performed. 
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Four  conditioning  indicators  were  used  to  evaluate  the 
stability  of  the  N  natrix.  These  were;  correlation  among 
the  paraneters,  snallness  of  the  determinant,  P-nuaier,  and 
the  overall  range  of  aagnitu'de  of  teras  in  the  N  matrix. 
Correlation  of  the  paraaeters  Mas  suggested  by  Ootila 
£  1975]-  The  idea  here  is»that  dependence  among  the  parame¬ 
ters  Mould  suggest  soae  dependence  in  the  columns  of  the  de¬ 
sign  aatrix.  The  saallness  of  the  determinant  and  the  P- 
nuaber  are  tests  used  by  Faddeev  and  Feddeeva  [15b3].  A 
singular  aatrix  has  a  deterainant  of  zero  and  can  not  te  in¬ 
verted,  but  a  aatrix  that  has  soae  dependence  can  be  invert¬ 
ed,  although  it  aay  be  ill-conditioned.  A  very  small  deter- 
ainant  suggests  that  the  aatrix  is  almost  singular.  Ibe 
P-nuaber  is  : 


|aax  \  j 

P-nuaber  =  - 

ain  \ 


vhere  /\  i  are  the  eigenvalues  of  the  H  aatrix.  Faddeev  and 
Faddeeva  suggest  that  a  very  large  P-ouaber  could  indicate 
ill-condition iag.  The  last  test  is  an  inspection  of  the  N 
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matrix  to  determine  if  there  is  a  large  range  in  the  magni¬ 
tudes  of  tens.  This  is  soaething  that  is  alaost  unavciau- 
ble  with  the  polynomial  aodel,  because  the  higher  order 
tens  of  the  polynoaial  are  the  products  of  variables  raised 
to  the  3rd  and  4th  powers.  The  best  way  this  author  has 
found  to  ainiaize  this  is  to  center  the  gndded  data  set. 
That  is,  put  the  origin  at  the  center  of  the  griu  system. 

The  "aatheaatical  sources"  of  error  are  truncation  and 
roundoff.  These  are  a  result  of  the  finite  representation 
of  nuabers  in  the  coaputer.  An  attempt  was  made,  in  this 
research,  to  avoid  or  ainiaize  these  errors  by  usn:^  aoutie 
precision  coaputation. 

The  "technical  source"  of  error  that  Ur  win  i  Ui.  WIN,  1^75] 
discusses  is  the  use  of  an  algoritha  that  is  poorly  written 
for  the  hardware  available.  This  should  not  it  a  ptoLiea 
because  a  proven  inversion  routine  (IHSL:  LINVil)  ana  a  very 
large  nainfraae  coaputer  (Aadahl  470V/6-11)  were  uses. 

Before  the  correlation  between  the  parameters  of  the  sur¬ 
faces  being  coapared  can  be  caculated,  the  parameters  must 
be  standarized.  This  is  because  there  is  a  large  range  in 
the  values  between  the  high  order  tens  and  the  low  cider 
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teras.  The  low  order  teras  are  large  values  cospared  to  the 
high  order  teras.  The  corresponding  order  teras  are  first 
grouped  then  standardized.  That  is,  the  ith  paraaeter  term 
of  each  trend  surface  is  standardized  in  relation  to  the  ith 
paraaeter  froa  all  other  trend  surfaces.  If  this  was  not 
done,  the  correlation  calculation  would  "amount  to  estimat¬ 
ing  the  variance  of  the  difference  corrected  for 

height.. . "[ NEBBIAH  and  SNEATH,  1966  ].  One  other  moaif icatiou 
of  the  paraaeters  is  necessary  before  the  correlation  is 
calculated,  the  oaission  of  the  first  paraeter.  This  is 
done  so  that  influence  of  the  constants  of  the  surface  (data 
ats)  will  not  be  present.  The  first  ter#  {2>  c  is  a  cons- 
nt  and  represents  the  intercept  of  the  surface  with  the  z 
axis.  The  second  ( Q  )  and  higher  k  J3n)  parameters  are  the 
ones  that  contain  the  surface  "character"  lururaaticn, 
therefore  these  are  the  paraaeters  that  are  ccm^ared. 

The  standardized  values  are  often  referred  to  as  "z" 
scores  in  statistics.  The  usual 
is  as  follows: 


x 


Ik 


X 


Ztk  * 


Sx 


aethod  of  standanzing  data 


x|K:  data  value  at  grid  point  jk. 

i  :  lean  value  of  *z*  values. 

s* :  the  standard  deviation  of  the  surface. 

A  final  consideration  is  the  fact  that  the  paraaeters  are 
not  independent  of  each  other.  This  can  he  seen  in  the  cor¬ 
relation  aatrix  of  the  paraaeters  figure  5.  The  high  order 
teras  of  the  polynoaial  are  correlated  with  the  first  order 
teras.  Because  the  correlation  coefficient  is  a  par aaete r ic 
statistic,  it  is  assuaed  that  saaples  (paraaeters)  are  inde¬ 
pendent.  The  independence  of  the  paraaeters  is  established 
through  the  use  of  orthogonal  polynoaials.  A  short  discus¬ 
sion  on  the  orthogonalization  process  follows. 

•••••cowuxATKHi  or  rutncTuamn 


The  correlation  aaong  pararateters  of  a  nonorthogonalizea 
bi-guartic  polynoaial  (25  teras). 
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figure  5:  Paraaeter  Correlation 
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The  parameter  estivates  tbat  are  generated  by  the  least 
square  solution  of  the  trend  surface  are  not  independent. 
Tbat  is  there  is  soae  correlation  between  the  parameter 
coefficients  in  the  design  aatriz  and  ultisitely  the  parame- 
ters  (estimates)  themselves.  Ibis  is  because  of  a  lack  of 
coaplete  independence  in  the  coluans  of  the  design  matrix, 
A,  for 


A 


=  Z 


the  general  aodel. 

A  standard  aetbod  used  to  transform  a  non-or thogonal  ma¬ 
trix  to  an  orthogonal  one  is  the  Graa-Schaidt  Process.  Cue 
can  think  of  the  coluans  of  A  as  a  basis  for  the  vector  of 
the  dependent  variables  Z.  By  applying  the  Gram-Schmiat 
Process  to  the  coluans  of  A,  the  independence  of  the  columns 
is  guaranteed.  For  example  let  <  a1  ,  a2  ,  a3,...,  a„  >  be  the 
coluans  of  the  design  aatriz  A.  Then  the  set  of  vectors  <  d1 
,  d2  ,  d3  , . . . ,  d„>  are  orthogonal  mben 


d1  =  a1  , 


Chapter  V 
iOTOCOBBBl ATI 01 


5.  1  IBTBOPOCTIOB 

"The  most  powerful  knowledge  of  a  surface  is  the  know¬ 
ledge  of  its  autocorrelation"  [ fiEUCKEB ;  1972  ]-  Autocorrela¬ 
tion  as  the  naae  suggests  is  the  correlation  of  something  to 
itself-  A  good  example  is  the  autocorrelation  of  a  wave- 
fora,  first  a  copy  of  the  the  nave  fora  is  superimposed  on 
itself.  Then  one  of  the  waveforms  is  shifted  a  small  dis¬ 
tance,  lag  (d) ,  and  the  waveforms  are  then  compared  along 
their  entire  coaion  length.  The  aoving  wavetoira  is  then 
shifted  again  so  that  the  total  displaceaent  is  2d,  the 
waveforas  are  then  coapared  again  and  the  process  is  contin¬ 
ued,  if  desired,  until  the  total  displacement  is  egual  to 
the  length  of  the  waveforas.  The  average  degree  of  similar¬ 
ity  between  all  points  a  certain  distance  apart  is  the  value 
that  is  calculated  with  autocorrelation.  The  autocorrela¬ 
tion  of  a  surface  can  be  used  as  a  representation  of  the 
surface.  By  coaparing  the  autocorrelation  functions  of  the 
gridded  data  sets,  the  original  data  sets  can  be  coapared. 
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Two  comparison  methods  that  use  the  autocorrelation  func¬ 
tion  aethod  are  applied  to  the  gridded  data  sets.  The  first 
method  uses  a  limited  autocorrelation  function,  in  which 
only  lags  along  specific  lines  are  used.  In  this  way,  the 
autocorrelation  along  cardinal  directions  is  the  only  con¬ 
sideration.  The  resultant  graphs  can  be  comparea  visually 
or  a  difference  curve  can  be  formed.  The  second  method  uses 
the  more  general  two-dimensional  spatial  autocorrelation 
function.  This  function  results  in  a  nap  or  gridded  data 
set.  By  comparing  the  autocorrelation  function  (nap)  of 
different  surfaces,  one  can  compare  the  onqmal  maps  or 
surfaces. 

The  theory  for  the  autocorrelation  function  comes  icon 
the  study  of  time  series.  Huch  of  the  terminology  reflects 
this  origin.  The  use  of  times  series  for  spatial  analysis 
is  well  developed,  and  there  are  many  good  examples  availa¬ 
ble  [  DAVIS;  1973, AGTEBBEHG;  1974  ].  The  logical  development 
for  spatial  autocorrelation  is  from  a  one  dimensional  analy¬ 
sis  to  the  multi-dimension  case.  The  theory  of  autocorela¬ 
tion  will  be  introduced  using  the  one  dimensional  case,  tol- 
lowing  that  a  discussion  of  spatial  autocorrelation  will  be 
presented. 
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5.2  §iuu±  an  Mjuagm  £±si 

The  following  discussion,  sections  5.2  and  5.3,  relies 
heavily  on  the  work  of  Agterberg  £1974].  Let  k=  -® 
«...  0,  1,... , o  for  soae  randoa  variable  x  of  an  ordered  ser¬ 
ies.  Aatocovariance  is  defined  as  : 

-  1  £  I  '  >  (  w  »  H  - 

fi  :  the  aean  of  the  ordered  series 
d  :  the  lag  or  displaceaent  along  the  series 

The  eapirical  autocovariance  function  is: 


The  autocorrelation  function  is: 


/Od  *  rd  /  ro  -  rd  /  6* 


the  variance  of  the  series 


4  1 


The  eapirical  case  is  [ 1GTEBBBBG; 1974 ]: 


r(d)  =  C(d)  /  C  (o) 


r  (d)  : 

autocorrelation  of 

lag 

(d) 

C(d): 

covariance 

at 

lag 

<<i) 

C  (o)  : 

covariance 

at 

lag 

zero 

(°) 

The  series  X  is  said  to  fce  "veakly"  stationary  when  it 
has  an  autocorrelation  function  that  is  constant  and  doesn't 
depend  on  location  along  the  series.  For  a  stationary  ser¬ 
ies 

r  — >  0  as  d  increases. 


I 


Bken  all  variables  X  have  the  saae  lean,  variance,  and  auto¬ 
correlation  function,  the  series  is  veak  stationary. 
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i  correlogram  is  a  sequence  of  autocorrelation  coeffi¬ 
cient st  r  with  d= 1 , 2,. . . , m.  The  figure  6  is  one  possible 
fora  of  a  correlograw. 


Figure  6:  Sample  Correlograw 


Before  woving  to  the  two  dimensional  case,  some  of  the 
important  properties  of  the  autocorrelation  function  will  be 
listed: 

1.  r(d)  <1  for  all  d 

2.  r(0)  =1 

3.  r  (d )  =  r(-d) 


The  autocorrelation  function  has  a  restricted  range  be¬ 
cause  it  uses  the  correlation  coefficient  to  compare  the 
series  as  it  is  displaced  by  different  lags.  Because  of  the 
restricted  range  of  the  function  it  seems  tc  be  a  good 
choice  for  map  comparison.  The  results  of  the  different 
methods  used  for  map  comparison  should  be  easier  to  ccmparc 
if  the  range  of  the  possible  results  are  the  same. 
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where  the  data  X  is  frow  an  (axn)  aatrix  of  regularly  gria- 
ded  data.  The  r  and  s  variables  describe  the  lag  direction 
and  distance.  For  exaaple,  C(2,2)  would  be  the  covariance 
of  all  points  that  are  seperated  by  2  grid  unite  over  (east) 
and  2  units  (south)  down.  The  first  eguation  will  compute 
the  covariance  of  all  points  in  the  directions  090°  -  180°  ; 
the  second  eguation  computes  the  covariance  of  alj.  lags  in 
the  180°-  270° directions.  C(0,0)  is  the  covariance  of  all 
points  at  zero  lag  distance. 
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The  autocorrelation  function  is: 


B  (r,s)  =  C(r,s)  /  C(0,0)  . 

This  result  can  be  used  to  fora  an  autocorrelation  aat. 
This  function  has  been  used  in  a  less  general  fora 
[LOON;  1976  ] 


C(d)  = 


n-d 

E 


<*l  >  <*.K  > 


ahere:  X,  :  the  noraalized  residuals  at  the  point 
n  :  the  nuaber  of  points 
d  :  the  lag  distance. 


This  is  the  saae  as  the  aore  general  function  if  the  condi¬ 
tion  B(r,0)  and  B(0,s)  Mere  u6ed  and  averaged  together. 
That  is,  only  the  M/S  and  E/B  lag  directions  are  considered. 
In  calculating  this  type  of  liaited  autocorrelation  func¬ 
tion,  one  could  be  lead  to  believe  that  the  function  was 
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rather  saooth  and  possibly  isotropic.  Plotting  average  pro¬ 
files  results  in  saooth  functions  that  are  siailar  to  figure 
7. 


Figure  7:  Liuited  Autocorrelation  Function 

Hhen  a  complete  autocorrelation  function  is  calculated 
and  averaged  to  give  a  profile,  the  result  is  rather  diffi¬ 
cult  to  interpret.  Figure  8  illustrates  an  autoccrrel ation 
function  for  a  complete  surface  evalutation,  without  regard 
for  spatial  orientation.  The  function,  calculated  in  its 
general  fora,  can  also  be  a  sap  where  r  and  s  au  the  coor¬ 
dinates.  That  is,  the  independent  variables  for  C(r,s)  and 
C(-r,s)  are  the  coordinates  for  the  autocorrelation  function 
values. 

Moraally,  the  calculation  is  not  carried  to  tne  largest 
possible  lag.  This  is  because  the  nuaber  of  cases  for  the 
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The  autococ relation  function  for  Borrison  Surface  III, 
calculated  without  considering  the  direction  of  each 
lag  distance. 


Figure  8:  Autocorrelation  without  respect  tc 
orientation 


larger  lags  is  not  sta tist ically  significant  and  the  corre¬ 
lation  value  can  exceed  one.  Figure  8  is  an  autocot  relation 
function  carried  to  the  largest  possible  lags;  it  was  calcu¬ 
lated  on  a  21x21  grid.  The  aatocorrelation  naps  for  April, 
flay,  and  June  have  correlation  isarithns  that  exceed  1.0. 
Because  it  is  possible  to  exceed  t.O,  care  should  be  taken 
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when  trying  to  interpret  the  autocorrelation  for  lags  that 
are  at  or  near  the  aap*s  aaxiaua  diaensions.  This  will  Le 
discussed  ia  chapter  eight  again. 


Chapter  FI 

GBIDDED  DATA  SUHFACES 


6.  1  IiTBODOCTIOM 


The  data  used  for  testing  and  evaluation  of  direct  corre¬ 
lation  and  trend  surface  prograas  case  form  two  artiticial 
surfaces  or  functions.  Once  the  programs,  used  for  map  com¬ 
parison,  were  considered  tested  and  operational,  pollution 
data,  provided  by  the  Ohio  Envicoaental  Protecti ox.  Agency 
(EPA) ,  was  evalutated  [Ohio  EPA,1979].  Both  types  or  data, 
artifical  and  real,  were  evaluated  in  a  21x21  gndded  for- 
aat.  Smaller  grid  sizes  were  possible  but  21x21  was  chosen 
because  it  provided  an  ovec-deterained  system  for  the  '.east 
square  fitting  of  the  polynoaials.  The  odd  size  grid  was 
chosen  to  aid  in  centering  the  data  sets.  Because  of  the 
large  nuaber  of  data  points  (441)  in  each  set  a  greater  de¬ 
gree  of  freedon  was  possible  when  calculating  the  correla¬ 
tion  coefficient  for  the  direct  correlation  aethod. 


6.2  UTIliat  SWiCIS 


Two  artifical  surfaces  were  used  for  the  initial  testing 
and  debugging  of  the  nap  comparison  prograas.  Both  qndaed 
data  sets  were  generated  by  functions  descrited  ty  Borasot 
[1971].  The  sinplest  surface  was  generated  by  the  exponen¬ 
tial  function 


2  =  exp[-<  X2  ♦  I1  )]- 


This  function,  called  Surface  11  by  Borrison,  was  used  to 
generate  data  values  for  a  21x21  grid.  The  z-values  were 
scaled  to  have  a  range  of  0  -  1000  units.  The  lasic  surface 
was  then  altered  by  adding  noise.  The  first  alteration  ad¬ 
ded  noise  signa-10,  the  second  added  signa-100,  and  the 
third  sigia=1000.  These  surfaces  thus  provided  gridueu  data 
that  was  correlated  by  varying  anounts.  The  generation  or 
noise  was  done  with  a  randon  nunber  generator  froa,  the 
Scientific  Subroutine  Package  (SSP) .  The  subroutine,  gauss, 
has  arguaents  that  require  a  nean  value  and  a  desired  stan¬ 
dard  deviation  for  the  randon  value  returned.  The  subrout¬ 
ine  returns  a  value  that  is  the  sun  of  the  aean  and  a  raado* 
value.  The  surfaces  were  altered  by  passing  each  of  the  or- 
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iginal  surface  data  points  to  GAOSS  as  a  mean  value.  GAUSS 
then  returned  a  nunber  that  was  the  sun  of  the  original  data 
value  and  a  randon  nunber  that  had  the  desired  standard  de¬ 
viation.  The  altered  surfaces  are  shovn  in  figure  5.  She 
original  surface  is  readily  visible  in  figure  9  t,  but  the 
original  surface  is  less  obvious  in  figure  9  c.  lhc  final 
alteration  with  signa=1000  caused  the  surface,  figure  s  d, 
to  be  unrecgonizable  with  little  of  the  original  trend  visi¬ 
ble. 


The  second  artifical  surface  was  generated  froir  a  tngo- 
nonetric  series  polynomial.  See  Horrison  [1971  ■  for  the  ex¬ 
act  equation.  This  complex  surface  is  referred  to  as  Sur¬ 
face  III  by  Horrison  [1971],  see  figure  10  a.  This  surface 
was  altered  also,  but  instead  of  adding  noise  it  was 
snoothed  using  a  nines  filter  [100N,1978].  The  nines  filter 
used  vas  based  on  the  Gaussian  Function,  ste 
[ LOOM ;  1978, p. 142  ]  for  a  coiplete  discription  of  the  filter 
used.  The  resulting  gridded  data  sets  are  similar  enough  to 
show  varying  correlation  between  then,  (see  figure  10b, c,d). 


li 
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Morrison  Surface  III  with  smoothings  from  a  nines  filter, 
a)  original  surface,  b)  first  smoothing,  c)  second 
smoothing,  d)  third  smoothing. 


Figure  10:  Horrisioa  Surface  III  with  smoothings. 
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6. 3  POLLOTIOl  SOBPACES 

There  are  a  number  of  pollutants  monitored  by  the  Ohio 
EPA,  Total  Suspended  Particulates  (TSP)  is  an  important  one. 
Substances,  both  solid  and  liquid,  that  are  emitted  into  the 
air  and  reaain  suspended  for  long  periods,  are  called  sus¬ 
pended  particulates.  Soae  of  the  coaion  forms  of  the  sus¬ 
pended  particulates  are  fly  ash,  process  dust,  fuses,  soot 
and  oil  aerosols  [Ohio  EPA,  1979  ].  The  reason  given  for  ao- 
nitoring  TSP  is  that  in  addition  to  Uniting  visifclity,  it 
also  daaages  the  human  respiratory  systea  ty  inter  f  err  mg 
with  the  lung's  natural  cleaning  process. 

The  saapling  of  TSP  is  accoaplished  by  drawing  large  vo¬ 
lumes  of  air  through  a  preweigbed  filter  to*,  a  specific 
time.  The  resultant  unit  is  aicrograms  per  cubic  liter  of 
air  passed  through  the  filter.  The  normal  sampling  duration 
is  a  twenty-four  hour  period  every  six  days.  The  sampling 
of  TSP  has  been  extensive  aud  continuous  since  the-  md 
1970’s,  thus  a  large  amount  of  data  was  readily  available. 

The  samples  used  to  construct  gridded  data  maps  of  ISP 
are  aonthly  averages  for  the  calendar  year  1990.  There  were 
a  total  of  60  continuously  operating  collection  sites  m  Cu¬ 
yahoga  and  surrounding  counties  for  1980.  Approxiaatly  half 
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of  the  sites  used  sere  located  in  the  City  of  Cleveland  with 
the  balance  of  the  sites  in  the  sarrounding  counties. 


The  box  uith  the  susbet  four  in  it  is  the  area  containing  the 
collection  stations,  it  corresponds  to  the  area  of 
the  nap  in  figure  12. 

Figure  11;  location  of  the  EEA  Monitoring  sites  used. 

The  obvious  clustering  of  the  data  presented  a  difficult 
problea  when  putting  the  data  in  a  21x21  grid  form.  The 
Geodetic  Science  Plotting  Package  (GSFF)  was  used  to  grid 
the  data.  GSPP  uses  a  weighted  average  of  the  nearest 
neighbors  to  predict  the  values  of  the  grid  points.  For 


this  pollution  data,  a  power  of  prediction  of  four  was  det- 
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erained  to  be  the  aost  suitable  value  to  use  for  the  weight¬ 
ing  function.  For  a  detail  discussion  of  the  weighted  aver¬ 
age  aethod  of  prediction  see  Sunfcel  [1980]. 

Using  the  aean  aonthly  TSP  values,  for  each  of  the  sixty 
sites,  twelve  gridded  data  aaps  were  constructed  using  the 
prediction  algoritba  of  GSPP.  See  appendix  &  fcr  ail  twelve 
aonthly  aaps.  The  isarithaic  aaps  produced  iron  tne  gridded 
data  represent  continuous  surfaces.  The  concentration  of 
TSP  is  the  unit  or  statistic  of  the  continuous  surface.  The 
drawing  of  the  isarithans  was  done  with  GSPP  and  figure  12 
is  an  exaaple  of  one  of  the  aonthly  aaps  produced. 

It  should  be  clear  that  the  aap  coaparison  techniques  used 
are  not  comparing  isarithaic  aaps,  but  rather  the  gridded 
data  sets  the  isarithaic  aaps  were  constructed  from. 


;is"c  3 

>:,3?e  -J 


46C9S 1  6  \  V. 

463, 53,  \  »c  ( 

460:C99  -  ,\M  <  _X 

4655330  -  A 

1459:566  ~Uk^  //  s 

45S55C5  -rSg^Vs  ®  '•'/  .  65 


45',s:6:  — • 


to 


'fm4h 


i3!  I  s*7 


I^CCIICB 

466.556  J 
455'53,  J 
4  5  5  5  ,  '  5  — | 


45,5546 
4  5  ?  r : 56 
4533,54 


4C0G06  41,59,  ,5656.'  445866  45'!$3  4’14?S  4eE'54 


Triangles  represented  collection  stations,  values  telow 
the  triangle  are  the  aonthly  averages  for  TSP 
at  the  sites  during  January. 


Figure  12;  Isarithaic  aap  for  January  I960  ISP. 


Chapter  VII 
11S0LTS 


Two  of  the  three  aethods  of  aap  comparison,  direct  cotre- 
iation  and  trend  surface,  were  applied  first  to  the  test 
surfaces  (figures  9,10)  and  finally  to  the  twelve  aonthiy 
TSP  surfaces  (see  appendix  A).  the  test  naps  were  expected 
to  be  correlated,  but  the  nature  of  the  pollution  surfaces 
was  unknown.  Comparison  of  the  sequential  months  was  loaned 
at  for  possible  groupings.  Ihe  relationships  tetween  all 
months  was  also  examined.  There  were  a  total  of  sixty-six 
possible  combinations  that  the  naps  could  bo  grouped  for 
comparison.  Spatial  autocorrelation  maps  were  termed  far 
each  of  the  twelve  aonthiy  TSP  maps.  Some  of  the  methods, 
notably  the  differences  of  the  limited  autocorrelation  func¬ 
tions,  were  generated  only  for  surfaces  that  had  shown  prem¬ 
ise  of  being  significantly  correlated.  A  large  amount  of 
data  and  maps  were  generated  by  the  different  comparison 
aethods.  For  a  complete  presentation  of  all  the  results  see 
the  appropriate  appendix.  Only  representative  examples  have 
been  included  in  this  discussion. 
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7.1  ABTIf IC11  SOEFACES 

The  artifical  surfaces,  figures  9  and  10,  were  evaluated 
by  applying  the  direct  methods  of  visual  comparison,  xsopach 
map  generation,  and  direct  correlation.  As  would  te  expect¬ 
ed,  the  artifical  surfaces  were  highly  correlated  visually. 
The  four  variations  of  the  Morrison  Surface  11  were  cos^-ared 
by  the  direct  correlation  method;  the  results  are  m  table 
1. 


r - 1 

TABLE  1  I 

Direct  Method  of  Correlation  for  Surface  II  1 

standard  deviation  of  noise  I 


0.0 

1.0 

100. 0 

1000.0 

o 

• 

o 

1.0 

— 

— 

— 

10.0 

.99 

1.0 

— 

— 

100.0 

.92 

.91 

1.0 

— 

1000.0 

.21 

.21 

.20 

1.0 

l 

It  can  be  seen  in  table  1  that  the  correlation  for  the  sur¬ 
faces  was  large.  The  correlation  between  the  smoothed  sur- 


faces  forced  fro*  Borrison  Surface  III  Here  also  high. 
There  were  known  variations  in  the  surfaces  but  the  correla¬ 
tions  indicated  little  variation  between  the  surfaces.  Ihe 
direct  correlation  aethod  didn't  appear  to  be  sensitive  to 
saall  changes  between  the  data  sets,  like  the  variations 
caused  by  adding  noise. 

The  isopach  aethod  was  applied  to  the  surfaces  ana  dif¬ 
ference  aaps  were  foraed.  The  result  was  a  set  of  that 

were  siailar  to  figure  1.  Figure  1  was  foraed  by  subtrac¬ 

tion  of  the  unaltered  Borrison  Surface  11  from  the  altered 
(Sigaa=100)  Surface  II.  The  result  was  an  ideal  random  sul- 
face,  without  any  trends  or  "topographic"  foras.  The  only 
difference  between  the  surfaces  was  noise  so  the  aaj.  (figure 
1)  has  a  raudoa  character.  Ihe  other  isopach  maps  (ora- 
had  a  siailar  appearence  with  different  isoline  values. 

Next,  the  trend  surface  aethod  was  applied  to  the  Morri¬ 
son  surfaces.  Ihe  correlation  between  the  parameters  or  the 
Surface  II  aaps  is  in  table  2.  Correlation  between  the  ori¬ 
ginal  surface  and  the  first  two  alterations  is  significant, 
but  there  is  a  negative  correlation  between  the  sigma- 1000 
and  the  other  surfaces.  This  relationhip  among  the  surface 
is  like  the  relationship  in  table  1.  The  sigaa=1000  surface 


is  such  less  correlation  with  the  rest  of  the  surfaces. 


This  is  because  the  range  of  the  unaltered  surfaces  was  1000 
and  the  noise  with  a  standard  deviation  of  1000  has  caused 
the  original  surface  to  be  radically  altered.  See  figure  H 
for  a  visual  verification  of  the  radical  change  in  the  sur¬ 
face. 


i - , 

I  TABLE  2  I 
J  Correlation  of  cubic  parameters  | 
|  Sigma  of  Noise  | 


o 

« 

o 

0.0 

1.0 

10.0 

100.0 

1000.0 

10.0 

.99 

1.0 

— 

-- 

100.0 

.79 

.76 

1.0 

— 

1000.0 

-.96 

-.95 

-.93 

1.0 

|  Correlation  between  parameters  of  a  cubic  polynomial  I 

I  fit  to  the  Horrison  Surface  II  and  alterations.  I 

I  1 

I  I 

i _ - _ _ _  , .  -  .  „  ... _ _ _ _ _ j 

Table  2  was  generated  by  comparing  the  parameters  of  a  cutic 
polynomial.  The  fit  of  the  polynomial  was  good  for  all  sur¬ 
faces,  see  table  3,  execpt  for  the  sigma=1000  surface. 
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I  I 

|  TABLE  3  | 

I  1 

|  Precentage  goodness  of  fit  foe  Surface  XI  I 


Sigaa  of  Surface  IX 

Goodness 

of  Fit,  % 

Cubic 

Bi-cubic 

Ei-guar tic 

0.0 

99.97 

100.0 

100.  0 

10.0 

99.78 

99.78 

99.80 

100.0 

83.  91 

84.14 

84.60 

1000.0 

6.8 

8.94 

16.60 

Table  4  is  the  result  of  using  a  bi-quartic  polynomial 
with  trend  surface  aethod.  Table  4  shows  a  drop  iu  correla¬ 
tion  between  the  sigaa=0,  sigaa=1Q,  and  sigma=10d  surf  acts. 
This  is  aost  likely  because  the  bi-quartic  polynomial  can  be 
aade  to  fit  the  gridded  data  sets  tetter  than  the  cubic  po- 
lynoaial.  The  better  fit  of  the  higher  order  polynomial 
causes  the  correlation  between  the  paraaeters  to  decease 
with,  an  increase  in  the  ability  of  the  polynomial  to  repre¬ 
sent  the  surface.  Hhat  seeas  to  be  happening  was  an  in¬ 
crease  in  sensitivity  of  the  trend  surface  method  corres¬ 
ponding  to  an  increase  in  the  ability  of  polynonial  aodel  to 
fit  the  gridded  data  sets.  A  siailar  trend  was  seen  with 
the  Borrison  Surface  ill  and  alterations.  The  size  ci  the 
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polynomial  used  should  improve  the  sensitivity  of  the 
method,  hut  the  overall  relationship  between  the  surfaces 
remains  the  same.  Pairs  of  surfaces  that  are  more  correlat¬ 
ed  will  remain  so  regardless  of  the  polynomial  model  used, 
this  can  be  seen  by  comparing  the  results  from  the  different 
models. 


TABLE  4 

Correlation  of  Bi-Quartic  Parameters 
Sigma  of  Noise 


O 

t 

O 

10.0 

100.0 

1000. 0 

o 

• 

o 

1.0 

— 

— 

— 

10.0 

CD 

CD 

• 

1-0 

— 

— 

100.0 

-79 

.69 

1.0 

— 

1000.0 

-.94 

-.90 

-.91 

1.0 

j 


Correlation  of  bi-guartic  polynomial  parameters. 


7-2  BUKI  flKIQJS  APPUSP  TO  I$f  lifS 


The  direct  visual  comparison  of  the  TSP  isanthmic  maps 
was  done  Hith  a  light  table.  The  similarity  between  each 
■ap  Has  determined  as  each  Has  overlayed  in  cosiination  with 
every  other  surface.  The  relation  between  the  temporally 
related  sooths  vas  examined  thoroughly  and  the  remaining  re¬ 
lationships  were  given  a  sore  superfical  evalution.  Ihe 
sixty-six  pairings  were  tedious  and  objectivity  was  ditti- 
cult  to  aaintain.  During  the  visual  comparison,  egual 
Height  Has  given  to  each  of  the  three  points  discussed  m 
chapter  3  under  the  heading  of  direct  visual  comparison.  An 
atteapt  Has  nade  to  construct  a  correlation  matrix  of  the 
pollution  naps,  but  the  results  were  not  consistent  or  cb- 
jective- 

The  direct  correlation  of  the  gridded  data  sets  without 
regard  for  spatial  location  was  accomplished  Ly  calculating 
the  correlation  coefficient  between  all  twelve  gndded  data 
sets.  This  resulted  in  a  synnetric  array  of  numbers  in 
which  the  diagonal  elenents  are  1.0,  a  perfect  correlation 
between  each  surface  and  itself.  The  off  diagonal  elements 
of  the  array  are  the  correlation  between  the  respective 
nonths.  Figure  13  is  the  correlation  between  the  monthly 
TSP  maps.  The  ten  highest  correlations  are  underlined.  As 
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can  be  seen  in  the  array,  there  are  aany  strong  to  moderate 
correlations  present.  Hone  of  the  surfaces  were  zero  or  ne- 
gativly  correlated.  The  correlation  appeared  to  be  guite 
strong  between  the  souths  that  are  temporally  associated. 
August  thrn  November  appeared  to  have  a  strong  correlation; 
they  were  all  related  with  correlations  greater  than  0.  d. 
The  strongest  correlations  were  Barch/Decemter  and  May/June. 
The  strong  to  noderate  correlation  throughtout  the  array 
suggests  that  all  TSP  surfaces  are  sinilar.  Host  likely  the 
lack  of  spatial  infornation  could  be  causing  everything  to 
seen  to  be  aore  correlated  than  it  actually  is. 

The  third  method  was  the  differencing  nethod  of  isopach 
naps.  The  TSP  maps  were  used  to  form  sixty-six  isopach 
maps,  of  which  fourteen  are  in  appendix  B.  The  fouteen  iso- 
pach  aaps  in  appendix  B  include  the  twelve  monthly  related 
isopach  aaps  and  three  other  aaps  that  will  be  discussed.  A 
notation  of  Jan/Feb  represents  the  isopach  map  formed  Ly 
subtracting  the  January  data  points  froa  the  Petuary  data 
points.  The  review  of  these  aaps  revealed  none  of  the  iso- 
pach  aaps  to  have  a  structure  or  appearance  similar  to  fig¬ 
ure  1.  Alaost  all  of  the  aaps  had  a  peak  or  depression  m 
the  upper  center  of  the  nap,  near  the  location  of  the  City 
of  Cleveland.  Ho  readily  apparent  nap  similarities  were 
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found  in  the  isopach  naps-  The  sources  of  the  correlation 
indicated  by  the  direct  nethod  were  not  apparent  xn  the  iso¬ 
pach  naps.  The  naps  that  had  strong  correlation  with  the 
direct  nethod  correlation  didn’t  appear  significantly  diffe¬ 
rent  than  the  naps  that  had  loner  correlations.  lhe  nap  of 
lowest  correlation,  Feb/Jun,  was  specifically  exanined  tut 
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it  had  no  special  features  either.  The  Apr/Jun  iscpach  nap 
did  sees  rather  flat,  which  sight  indicate  that  the  surfaces 
of  April  and  June  Here  similar,  but  there  was  a  strong  de¬ 
pression  near  Cleveland  on  the  sap,  suggesting  there  were, 
also  some  sigiuificant  differences  between  Bay  and  June. 

7.3  TREND  SURFACE  ajTHQD  APPLIED  10  1J3JP  HAPS 

The  monthly  gridded  data  sets  ot  TSP  were  fitted  with 
three  different  polynomial  models.  A  cubic,  n-cubic,  a:.d  a 
bi-g  uaLtic  polynomial  were  used  as  mathematical  monels.  Ihe 
design  matrix  in  each  case  was  or thogonalized  sc  that  the 
estimates  of  the  parameters  would  be  independent.  Ihe  (.di¬ 
ameters  were  then  standardized  by  the  method  described  in 
chapter  4.  The  resultant  standardized  parameters  rci  each 
surface  were  coapareu  using  the  correlation  cot f Licient.  An 
array  was  formed  with  the  resultant  correlations ,  (see  fig¬ 
ure  1 4 )  . 

The  goodness  of  fit  for  the  higher  order  polynomials  was 
■  uch  better  than  the  lower  order  polynomials ,  set  tails  ‘i. 
The  sum  of  squares  tor  the  residuals  <V  V>,  deettased  also 
with  the  larqer  polynomials.  This  would  seem  to  indicate 
that  the  mathematical  model  (polynomial)  was  tecomiug  flexi¬ 
ble  enough  to  fit  the  irregular  TSP  surfaces.  If  the  sun- 
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TABLE  5 

Percentage  Goodmess  of  Fit  for  TSP  Baps 


Goodness  of  fit,  * 


Monthly  TSP  Bap 

Cubic 

Bi-cutic 

Bi-guartic 

Jan 

40.21 

45.00 

58.52 

Feb 

62.98 

68.29 

75.25 

Bar 

28.68 

32.90 

52.43 

Apr 

27.90 

29.56 

46.34 

flay 

27.09 

28.78 

45.75 

Jun 

28.45 

30-26 

48.72 

Jul 

31.60 

39.56 

57.  43 

Aug 

31.64 

34.85 

54.  4  1 

Sep 

35.27 

39.26 

54.  1  1 

Oct 

26.26 

30.37 

47.95 

Nov 

39.35 

46.28 

55.21 

Dec 

39.53 

44.92 

bl.CI 

faces  were  s aoother  the  polyaoaial  aodels  would  have  fitted 
even  better. 

The  correlation  Matrix  foraed  gives  an  indication  of  the 
relative  similarity  between  all  the  compared  surfaces.  The 
surfaces  that  had  strong  correlation  with  the  direct  method 
showed  less  correlation  when  the  trend  surface  method  was 
applied.  The  correlations  of  figure  13  are  larger  than 
those  of  the  the  trend  snrface  aethods  see  figure  14  and  ay- 
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Correlation  of  TSP  aaps  using  the  trend  surface  aetnod  with 
a  bi-guartic  polynomial. 

Figure  14:  Correlation  natrix  for  bi-guartic  model 

pendix  C.  There  are  negative  and  near  zero  correlations 
where  large  correlations  had  appeared  with  the  direct  meth- 
od.  The  range  of  correlation  values  has  increased  signifi¬ 
cantly  with  the  trend  surface  aethod.  Ihe  two  highest  value 
pairings,  flay/June  and  Barch/Dcceaber,  renained  the  hignest 


values 


b'i 


The  shift  of  values  between  the  trend  Methods  of  diffe¬ 
rent  aatheaatical  Models,  aost  likely  was  a  result  cl  the 
iaproved  fitting  of  the  polyooaial  to  the  data.  That  is, 
with  iaproved  fitting  of  the  aodel,  the  coaparison  ol  esti- 
aated  paraaeters  would  give  a  wore  accurate  comparison  of 
the  TSP  gridded  data  aaps.  See  Appendic  C  for  the  correla¬ 
tion  of  the  TSP  sufaces  using  a  cubic  and  bi-culic  mathemat¬ 
ical  models.  The  correlation  array  for  the  biyuartic  model 
should  be  the  aost  accurate  relation  of  the  three  models 
used,  because  the  bi-quartic  aodel  resulted  in  the  Lest 
goodness  of  fit  and  lowest  <v,v>  value. 

7. H  hOTOCOMlghlTIOl  BKtHODS 

The  autocorrelation  function  was  used  to  compare  gridded 
data  aaps  in  two  ways.  The  first  aethod  used  a  limited  au¬ 
tocorrelation  function,  by  which  the  general  spatial  auto¬ 
correlation  function  (see  section  5.3)  was  evaluted  in  only 
one  direction.  Autocorrelation  in  the  North/South  direction 
was  used  to  fora  a  series  of  curves  with  B(0,s)  form  of  the 
general  autocorrealtion  function.  An  B(0,  s)  function  was 
calculated  for  each  of  the  twelve  TSP  aonthly  maps.  Then 
sizty-six  difference  curves  were  foraed,  with  print ei  plots, 
froa  the  aonthly  curves.  The  printer  plots  were  checked  to 
see  if  the  surfaces  that  had  showed  correlation  Ly  a  previ- 
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ous  aethods  would  appear  different  froa  the  rest  ot  the 
curves.  The  (larch  and  Deceaber  carves  were  similar  and  re¬ 
sulted  in  a  difference  curve  that  was  nearly  zero.  Ihe  Hay 
and  June  curves  were  not  as  siailar  as  one  would  expect 
based  on  the  results  of  the  other  aethods.  looking  at  the 
plots  of  the  aonthly  curves,  in  pairs,  seemed  sore  useful 
than  looking  a  the  curve  differences. 

The  aajor  problea  with  the  liaited  autocorrelation  method 
was  that  by  choosing  a  particular  profile  to  work  with,  the 
rest  of  the  autocorrelation  inforaation  was  neglected. 

Piguce  15,  shows  the  profile  for  the  Sorth/South,  hast/ West, 
Sorthwest/Southeast,  and  Hortbeast/Southwest  directions  on 
the  autocorrelation  aap  of  January.  They  are  differ  tut 
enough  that  to  use  one  as  a  representation  of  the  surface, 
auch  additonal  autocorrelation  inforaation  would  te  ignored. 
Possibly,  a  potential  correlation  would  be  aissed  this  way. 

The  second  aethod  of  coaparison  used  the  general  fora;  of 
the  function.  The  twelve  ISP  gridded  data  sets  were  evalu¬ 
ated  and  the  autocorrelation  aap  of  each  ISP  surface  was 
foraed.  The  general  autocorrelation  function  would  produce 
oaly  one  half  of  the  autocorrelation  aap,  but  that  was  all 
that  was  necessary  because  of  the  syaaetry  of  the  function. 
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p. ~\  f?c 


Figure  15:  Autocorrelation  profiles  for  January 

la  order  to  fora  an  autocorrelation  sap, (see  figure  ibj, 
aatrices  that  were  generated  by  the  general  spatial  autocor¬ 
relation  were  put  together  to  fora  the  lover  half  c t  the 
aap.  The  coiaon  boundary  between  these  aatrrccs  caused  the 
eliaenation  of  one  coloan,  thus  the  two  (r  x  s)  arrays 
foraed  a  single  (s  x  2 r-1)  array.  That  was  only  one  bait  of 
total  autocorrelation  function  aap.  The  second  half  was 
foraed  by  rotation  of  the  first  half  by  180  degrees  aLcut 
the  east/west  axis  and  180  degrees  about  the  north/south 


Figure  it>:  January  autocorrelation  aar. 


axis.  There  Mill  be  a  conon  con  shared  by  the  upper  and 
lower  halves  when  joined,  thus  causing  a  loss  cf  one  row  in 
the  nap-  The  final  nap  array  is  (2s) -1  by  (2r)-l  ,  with  the 
northwest  guadrant  a  synnetric  inage  of  the  southeast  quad¬ 
rant.  1  siailar  syaaetry  was  true  between  the  southwest  and 
northeast  quadrants.  The  original  TSP  naps  were  (21  x  21) 
the  resultant  autocorrealtion  naps  were  (41  x  41). 
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The  syaaetry  of  the  aap  is  evident  in  figure  16.  The 
center  of  the  aap,  point  (21,21),  is  where  a  correlation  of 
1.0  should  exist.  The  direction  and  distance  out  from  the 
center  yields  the  autocorrelation  for  a  siailar  shift  and 
direction  on  the  TSP  surface.  The  boundary  areas  of  the  au- 
tocorrealtion  aap  are  subject  to  large  variations.  At 
tiaes,  the  correlation  value  will  exceed  1.0  or  -1.0,  (see 
figure  16).  This  is  aost  likely  caused  by  the  suulei  of 
saaples  available  for  lag  distances  near  the  extremes  of  the 
surface.  The  saaller  lags  have  lore  saaples  available,  thus 
the  saaple  is  a  aore  accurate  representation  of  the  popula¬ 
tion  it  is  intended  to  represent. 

Because  of  the  large  variation  near  the  edge  or  the  auto¬ 
correlation  aap,  the  values  aay  not  be  as  reliable  as  the 
ones  near  the  center  of  the  aap.  One  should  use  the  center 
2/3  of  the  diaension  of  the  aap.  Figure  8  doesn't  have 
correlation  values  exceeding  1.0  if  the  lags  less  than  2/3 
of  the  surface  are  used.  Inspection  of  the  spatial  autocor¬ 
relation  aaps  (appendix  0)  reveal  that  none  of  the  surfaces 
have  isolines  (correlation)  greater  than  1.0  if  only  the  in¬ 
ter  2/3  of  the  aap  is  considered.  For  aap  comparison,  only 
the  center  part  of  the  function  should  be  used.  By  avoiding 
the  large  boundary  variations,  a  aore  accurate  comparison 


should  be  possible.  To  cospare  these  surfaces  is  somewhat 
easier  than  cosparing  the  original  TSP  gridded  data  saps  be¬ 
cause  the  autocorrealtion  surfaces  are  symmetric.  With  a 
large  nuaber  of  surfaces  to  cospare,  the  problem  ma  y  still 
be  present  that  the  investigator  say  have  to  judge  between 
sany  possible  sap  cosbinations. 

The  fora  of  the  autocorrelation  function  was  similar  for 
each  aap,  (see  appendix  C  for  all  the  autocorrelation  maps). 
The  center  area  is  peaked,  with  the  value  at  or  near  1.0  and 
a  symmetric  cone  dropping  away  on  all  sides  to  a  value  or  . 5 
to  -2-  A  completely  isotropic  surface  would  be  plotted  as 
concentric  contours.  The  isotropic  nature  of  the  T3F  sur¬ 
face  was  lost  after  4-5  lags  from  the  center  in  any  direc¬ 
tion.  This  would  seen  to  be  a  strong  indication  that  the 
gridded  pollution  data  used  in  this  study  is  not  isotropic. 

A  visual  comparison  of  the  two  pairs  of  strongly  corre¬ 
lated  nonths,  Hay/June  and  Barch/Eecember ,  revealed  that  the 
autocorrelation  naps  were  also  sinilar.  It  was  also  found 
that  Hay/Jume/July  have  very  sinilar  autocorrelation  naps. 
To  avoid  the  problea  of  comparing  many  different  combina¬ 
tions  of  naps,  the  trend  surface  method  were  applied  tc  the 
autocorrelation  maps.  A  bi-guartic  polynomial  was  applied 


Chapter  V2I1 
COICLOSIOI 


Three  general  methods  of  sap  comparison  have  keen  pre¬ 
sented;  direct  methods,  trend  surface  and  autocorrelation. 
Each  aethod  compares  the  gridded  data  surface  in  a  different 
aanner  and  thus  will  give  a  different  indication  of  similar¬ 
ity  of  the  naps  compared.  The  relation  between  saps,  if 
looked  at  with  only  one  aethod,  aay  not  give  the  investiga¬ 
tor  a  complete  evaluation  of  the  similarity  between  surfac¬ 
es.  Bore  than  one  aethod  should  be  used  to  evaluate  the  si¬ 
milarity  between  surfaces.  To  investigate  a  large  numLer  of 
maps,  a  testing  strategy  should  be  used. 

An  initially  guick  and  inexpensive  aethod  cf  map  compari¬ 
son  is  the  test  of  correlation  without  regard  to  spatial  lo¬ 
cation.  The  aap  similarity  will  be  somewhat  exaggerated 
with  this  test,  because  the  surfaces  will  appear  more  corre¬ 
lated  than  they  really  are.  This  happened  with  the  test 
data  and  the  aonthly  TSP  surfaces.  They  all  showed  high 
correlations.  If  sore  investigation  is  warranted,  theD  the 
next  test  could  be  the  trend  surface  aethod.  If  the  direct 
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correlation  aethod  reveals  very  high  correlation  between  the 
aaps,  then  the  direct  aethod  aay  not  te  sensitive  enough  to 
sort  oat  the  differences  in  the  set  of  saps  the  investigator 
is  using.  If  the  investigator  feels  that  the  comparison  us¬ 
ing  the  direct  aethod  is  not  seperating  the  nap  set  suffi¬ 
ciently  the  trend  surface  aethod  should  be  the  next  aethod 
applied. 

For  this  aethod,  a  large  bi-polynoaial  should  be  used  if 
the  surface  is  highly  irregular.  In  the  original  use  of 
this  aethod  by  Rerriaa  and  Sneath  [1966]  a  cubic  polynomial 
was  the  largest  polynoaial  used  because  the  surfaces  could 
be  fitted  well  without  going  to  larger  polynomials.  This 
investigation  required  a  aucb  larger  polynoaial  to  get  a 
reasonable  goodness  of  fit,  see  table  5.  Once  a  polynomial 
aodel  is  chosen,  then  the  surface  can  be  tested.  An  orthog¬ 
onal  polynoaial  should  be  used  to  aaintam  independence 
aaong  the  estiaated  paraaeters  that  are  to  be  tested  tor 
correlation,  see  figure  5.  Correlation  testing  kith  trend 
surfaces  aethods  reguires  considerably  aore  numerical  pro¬ 
cessing  than  the  direct  aethod,  but  the  resultant  correla¬ 
tion  should  give  a  good  indication  of  the  siailarity  in 
shape  between  the  surfaces  that  are  coapared. 
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At  this  point  in  the  testing  sequence,  there  nay  be 
strong  correlations  evident  from  either  the  direct  nethcd  or 
the  trend  surface  Method.  Bj  choosing  the  aost  promising  or 
interresting  relationship  froa  the  correlation  latia  and 
then  applying  the  visual  aethod  the  investigator  would  avoid 
studying  every  single  nap  comparison-  The  iscpach  map  of 
the  surfaces,  that  have  an  indication  of  strong  correlation 
with  direct  or  trend  surface  methods,  could  also  be  formed 
and  looted  at,  thus  avoiding  the  construction  of  many  iso- 
pach  aaps. 

If  the  investigator  wished  a  nore  thorough  investigaticn, 
the  autocorrelation  of  the  surfaces  should  be  the  next  step 
to  study  the  siailarity  of  the  aaps  being  compared.  The 
calculation  of  the  liaited  autocorrelation  function  seems  to 
have  a  minimal  value,  because  of  its  inablity  to  represent 
non-isotropic  surfaces,  see  figures  13  and  15.  The  general 
autocorrelation  function  aap  seeas  to  be  good  method  of 
representing  a  surface  in  an  unigue  to  facilate  the  visual 
comparison  of  gridded  data  maps.  The  calculation  of  the 
general  spatial  autocorrelation  function  aap  requires  more 
nuaerical  effort  than  the  tread  surface  method,  but  the  ef¬ 
forts  seen  worthwhile.  The  quickest  way  to  evaluate  the 
correlation  between  aany  autocorrelation  aaps  is  to  use  the 
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tread  surface  aethod  applied  to  the  autocorrelation  saps. 
Ac  advantage  in  comparing  the  autocorrelation  naps,  over 
coaparing  the  original  data,  is  the  syaaetry  of  the  correla¬ 
tion  aap.  Additionally,  the  better  goodness  of  fit  possible 
with  the  syaaetric  surface  improves  the  accuracy  of  the  cor¬ 
relation  coaparison.  Visual  coaparison  is  also  improved  te- 
cause  of  the  syaaetry  of  the  aaps,  see  figure  16.  The  use 
of  visual  coaparison  is  best  reserved  for  the  final  check  of 
aaps  that  are  to  be  coapared.  Unless  the  number  of  maps  to 
be  coapared  is  saall  the  visual  aethod  should  te  used  last. 

To  coapare  aaps  by  only  one  of  the  techniques  presented 
here,  aay  cause  the  investigator  to  overlook  important  simi¬ 
larities  or  differences  that  are  present.  The  Lest  strategy 
for  aap  coaparison  using  the  aethods  outlined  in  this  paper, 
appears  to  be  the  use  of  a  battery  of  tests  that  have  in¬ 
creasing  coaputational  overhead,  that  the  investigator  can 
follow  as  far  as  is  necessary.  It  is  felt  Ly  the  author 
that  no  one  test  is  going  to  tell  the  full  story  about  the 
siailarity  between  a  set  of  gidded  data  naps,  but  in  concert 
with  each  other  the  aethods  can  be  adeguateiy  used  to  ana¬ 
lyze  a  set  of  gridded  data  naps. 


Appendix  A 

ISABITBBIC  MAPS  OP  TBE  HOITBLT  ISP  DAI A 


This  appendix  contains  the  isarithaic  representation  of 
the  gridded  data  naps  of  Total  Suspended  Particulates.  2 he 
TSP  naps  are  of  the  aonthly  averages  of  the  pollutants  from 
sixty  aonitoring  sites  in  Borthern  Ohio-  2he  isarithaic 
aaps  are  foraed  on  a  21  x  21  grid. 
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IppMdii  C 

TBSIO  S0BF1CE  P1BBBBTEB  COBBBL1TIOB  BESOLTS 


The  paraaeters  of  different  size  polynoaial  were  used  for 
coaparison  of  the  TSP  nap.  k  cubic  (9  ten)  ,  li-cufcic(16 
tern),  aad  bi-guartic  (25  tera)  polynoaial  were  each  usea  as 
aatheaatical  aodels  for  the  trend  surface  aethcd.  The  re¬ 
sults  of  fitting  a  cubic  and  bi-cubic  polynoaial  to  the  15? 
aaps  and  the  calculation  of  the  correlation  between  the  par- 
aaetrs  of  each  of  the  12  ISP  aaps  are  in  this  appendix.  1  he 
result  of  each  coaparison  is  a  12x12  syaaetric  array  with 
correlation  between  each  possible  pairing  of  the  ISP  monthly 
pollution  aaps.  The  result  of  the  bi-guartic  polynomial  are 
in  figure  14.  The  third  table  in  this  appendix  is  the  cor¬ 
relation  between  the  autocorrelation  aaps  that  were  formed 
froa  the  twelve  TSP  aaps.  The  bi-guartic  polynomial  was 
used  to  coapare  the  spatial  autocorrelation  aaps. 
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Appendix  0 
AOTOCOBBU.ATIOB  BAPS 


This  appendix  contains  the  autocorrelation  naps  for  the 
twelve  aonthly  TSP  surfaces.  For  each  nonthly  1st  nap  there 
is  one  isarithaic  aap  and  one  surface  representaion.  Cn  the 
isarithaic  aaps  the  syaaetry  of  the  autocorrelation  is  read¬ 
ily  apperant.  The  B/8  gua-l^ant  is  syaaetnc  with  respect  to 
the  S/E  quadrant.  These  aaps  are  produced  cn  a  41  x  41 
grid.  The  saall  triangles  ace  grid  intersections  where  data 
values  were  predicted. 
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TWO  -  DIFIENSIONFIL  RUTOCORREIRTIQN  FUNCTION 
.FOR*  FEBBO  TSP  GRIDDED  DRTR  ............ 

ISRRITHBS  IN  CORRELRTIDN  UNITS  *  M1XM1  • 


DIRECTION  or  UIEU 
LONGITUDE  -  -60’  00 ' 

LRTITUOE  -  30'  00' 
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,TLJO-  DIRENSIONHL  flUTOCORRELRTION  FUNCTION 

FDR*  RPRSO  TSP  GRIDDED  DRTR  . . 

lTSRRITHPIS  IN  CDRRELRT ION  UNITS  »  41X41  . 
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TUO-DinENSIONflL  AUTOCORRELATION  FUNCTION 

FOR-  P1RY80  TSP  CRIDDED  DRTR  •  ••■ . 

ISRRITHRS  IN  CORRELATION  UNITS  *  M1X41  - 


jTUO-DinENSIONRL  AUTOCORRELATION  FUNCTION 

FOR*  JUN80  TSP  CRIOOEO  DRTR  . . 

iISRRITHRS  IN  CORRElRTION  UNITS  ■  M1XM1  » 
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,TUC-OI'r.ENSICNRL  RUTOCORRElRTION  FUNCTION 
FOR-  JULBO  TSP  GRIOOED  OR TR  ............ 

ISRRITHFIS  IN  CORRELATION  UNl T 5  *  4 1X41  * 
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rrwo-DinENsiONRL  rutocorrelrtion  function 

iFOR*  RUG80  TSP  GRIQDED  DRTR  . 

ISRRITrlR5  IN  CORRELATION  UNITS  *  41X41  « 
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TUO-DIFIENSIONflL  AUTOCORRELATION  FUNCTION 
FOR-  OCTBO  TSP  GRIOOED  DRTR  ............ 

ISRRITHRS  IN  CORRELRT ION  UNITS  ■  HIXMl  * 
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